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Abstract 

An electron spin moving adiabatically in a strong, spatially non-uniform 
magnetic field accumulates a geometric phase or Berry phase, which might 
be observable as a conductance oscillation in a mesoscopic ring. Two contra- 
dicting theories exist for how strong the magnetic field should be to ensure 
adiabaticity if the motion is diffusive. To resolve this controversy, we study 
the effect of a non-uniform magnetic field on the spin polarization and on the 
weak- localization effect. The diffusion equation for the Cooperon is solved ex- 
actly. Adiabaticity requires that the spin-precession time is short compared 
to the elastic scattering time — it is not sufficient that it is short compared to 
the diffusion time around the ring. This strong condition severely complicates 
the experimental observation. 

PACS: 03.65.Bz, 72.80.Ng, 73.20.Fz, 73.23.-b 



Typeset using REVTeX 



1 



I. INTRODUCTION 



The adiabatic theorem of quantum mechanics implies that the final state of a particle 
that moves slowly along a closed path is identical to the initial eigenstate — up to a phase 
factor. The Berry phase is a time-independent contribution to this phase, depending only 
on the geometry of the path.0 A simple example is a spin-1/2 in a rotating magnetic field B, 
where the Berry phase equals half the solid angle swept by B. It was proposed by Stern! to 
measure the Berry phase in the conductance G of a mesoscopic ring in a spatially rotating 
magnetic field. Oscillations of G as a function of the swept solid angle were predicted, similar 
to the Aharonov-Bohm oscillations as a function of the enclosed flux.i 

An important practical difference between the two effects is that the Aharonov-Bohm 
oscillations exist at arbitrarily small magnetic fields, whereas for the oscillations due to the 
Berry phase the magnetic field should be sufficiently strong to allow the spin to adiabatically 
follow the changing direction. Generally speaking, adiabaticity requires that the precession 
frequency uj-q is large compared to the reciprocal of the characteristic timescale t c on which 
B changes direction. We know that uj-q = gfi^B/2h, with g the Lande-factor and /zb the 
Bohr magneton. The question is, what is t c l In a ballistic ring there is only one candidate, 
the circumference L of the ring divided by the Fermi velocity v. (For simplicity we assume 
that L is also the scale on which the field direction changes.) In a diffusive ring there are two 
candidates: the elastic scattering time r and the diffusion time around the ring. They 
differ by a factor r^/r ~ (L/£) 2 , where I = v r is the mean free path. Since, by definition, 
L ^> i in a diffusive system, the two time scales are far apart. Which of the two time scales 
is the relevant one is still under debate.@ 

Stern's original proposal! was that 

w B »- (1.1) 

T 

is necessary to observe the Berry-phase oscillations. For realistic values of g this requires 
magnetic fields in the quantum Hall regime, outside the range of validity of the semiclassical 
theory. We call Eq. ( LTD the "pessimistic criterion" . In a later work,! Loss, Schoeller, and 



Goldbart (LSG) concluded that adiabaticity is reached already at much weaker magnetic 
fields, when 




(1.2) 

This "optimistic criterion" has motivated experimentalists to search for the Berry-phase 
oscillations in disordered conductors,^ and was invoked in a recent study of the conductivity 
of mesoscopic ferromagnets.i In this paper, we re-examine the semiclassical theory of LSG 
to resolve the controversy. 

The Berry-phase oscillations in the conductance result from a periodic modulation of 
the weak-localization correction, and require the solution of a diffusion equation for the 
Cooperon propagator. To solve this problem we need to consider the coupled dynamics of 
four spin-degrees of freedom. (The Cooperon has four spin indices.) To gain insight we first 
examine in Sec. the simpler problem of the dynamics of a single spin variable, by studying 
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the randomization of a spin-polarized electron gas by a non-uniform magnetic field. We 
start at the level of the Boltzmann equation and then make the diffusion approximation. 
We show how the diffusion equation can be solved exactly for the first two moments of the 
polarization. The same procedure is used in Sec. [Ill] to arrive at a diffusion equation for 
the Cooperon. This equation coincides with the equation derived by LSG in the weak-field 
regime oj^t <C 1, but is different in the strong-field regime oj^t > 1. We present an exact 
solution for the weak-localization correction and compare with the findings of LSG. 

Our conclusion both for the polarization and for the weak-localization correction is that 
adiabaticity requires wbt 3> 1- Regrettably, the pessimistic criterion ( |1.1| ) is correct, in 
agreement with Stern's original conclusion. The optimistic criterion ( |1.2| ) advocated by 
LSG turns out to be the criterion for maximal randomization of the spin by the magnetic 
field, and not the criterion for adiabaticity. 

II. SPIN-RESOLVED TRANSMISSION 

A. Formulation of the problem 

Consider a conductor in a magnetic field B, containing a disordered segment (length L, 
mean free path I at Fermi velocity v ) in which the magnetic field changes its direction. An 
electron at the Fermi level with spin up (relative to the local magnetic field) is injected at 
one end and reaches the other end. What is the probability that its spin is up? 



B 




FIG. 1. Schematic drawing of a two-dimensional electron gas in the spatially rotating magnetic 
field of Eq. Q, with / = 1. 

For simplicity we take for the conductor a two-dimensional electron gas (in the x-y plane, 
with the disordered region between x = and x = L), and we ignore the curvature of the 
electron trajectories by the Lorentz force. The problem becomes effectively one-dimensional 
by assuming that B depends on x only. We choose a rotation of B in the x-y-pl&ne, according 
to 
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B 



(x, y, z = 0) = (B sin 77 cos ^f-, B sin 77 sin S cos ; 



(2.1) 



with t] and / arbitrary parameters. The geometry is sketched in Fig. |I|. We treat the orbital 
motion semiclassically, within the framework of the Boltzmann equation. (This is justified if 
the Fermi wavelength is much smaller than £.) The spin dynamics requires a fully quantum 
mechanical treatment. We assume that the Zeeman energy g^B is much smaller than the 
Fermi energy |mt> 2 , so that the orbital motion is independent of the spin. 

We introduce the probability density P(x,<f),£,t) for the electron to be at time t at 
position x with velocity v = {y coscp, v sin 0, 0), in the spin state with spinor £ = (£1,^2)- 
The dynamics of £ depends on the local magnetic field according to 



dt 



2h 



(2.2) 



where cr = (a x , <j y , a z ) is the vector of Pauli matrices. It is convenient to decompose £ 
Xift + X2^i into the local eigenstates of B • cr, 



cos f e~ i7Tfx/L 
sin f e infx/L 



- sin \ e -™f x / L 
cos f e iwfx/L 



(2.3a) 



B- 0-^ = 5^, B-<t( 



(2.3b) 



and use the real and imaginary parts of the coefficients X15X2 as variables in the 
Boltzmann equation. The dynamics of the vector of coefficients c = (01,02,03,04) = 
(Rexi,Imxi,Rex2,Imx 2 ) is given by 



dc 



T 



-Mc, M = Mq + Mi cos ( 



(2.4a) 
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(2.4b) 



where uj-q = gfisB /2h is the precession frequency of the spin. The Boltzmann equation takes 
the form 



d . 



cos. 



dP 

dx 



(2.5) 



where we have assumed isotropic scattering (rate 1/r = v/t). 

We look for a stationary solution to the Boltzmann equation, so the left-hand-side of Eq. 
( |2.5| ) is zero and we omit the argument t of P. A stationary flux of particles with an isotropic 
velocity distribution is injected at x — 0, their spins all aligned with the local magnetic field 
(so X2 = at x = 0). Without loss of generality we may assume that %i = 1 at x = 0. No 
particles are incident from the other end, at x = L. Thus the boundary conditions are 
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P(x = 0, 0, c) = 8(ct - l)5(c 2 )5(c 3 )5(c 4 ) if cos > 0, 
P(x = L, 0, c) = if cos < 0. 



(2.6a) 
(2.6b) 



This completes the formulation of the problem. We compare two methods of solution. 
The first is an exact numerical solution of the Boltzmann equation using the Monte Carlo 
method. The second is an approximate analytical solution using the diffusion approximation, 
valid for L ^> £. We begin with the latter. 



B. Diffusion approximation 

The diffusion approximation amounts to the assumption that P has a simple cosine- 
dependence on 0, 

P(x, 0, c) = jV(ar, c) + J(x, c) cos 0. (2.7) 

To determine the density N and current J we substitute Eq. (|2.7|) into Eq. (|2.5|) and integrate 
over 0. This gives 

81 8 

e— = -—(2M cN + M lC J). (2.8) 
ox oc 

Similarly, multiplication with cos0 before integration gives 

dN d 

£ ~dx~ = ~d~c (M ° cJ + MlCiV) " 1 (2 ' 9) 

Thus we have a closed set of partial differential equations for the unknown functions N(x, c) 
and J(x,c). Boundary conditions are obtained by multiplying Eq. (|2.6| ) with cos0 and 
integrating over 0: 

N(x = 0, c) + | J(x = 0, c) = 8(d - l)5(c 2 )5(c 3 )5(c 4 ), (2.10a) 
N(x = L,c)- jJ(x = L,c) = 0. (2.10b) 

We seek the spin polarization p = cf + c\ — c\ — c| of the transmitted electrons, charac- 
terized by the distribution 

_ Jdc J{x = L, c)5(c\ + c 2 2 -c 2 3 -cl-p) 
[P) ~ JdcJ(x = L,c) • 1 j 

(The notation Jdc = Jdc\ Jdc^ Jdc^ Jdc^ indicates an integration over the spin variables.) 
We compute the first two moments of P(p). The first moment p is the fraction of transmitted 
electrons with spin up minus the fraction with spin down, averaged quantum mechanically 
over the spin state and statistically over the disorder. The variance Var p = p 2 — p 2 gives 
an indication of the magnitude of the statistical fluctuations. 

Integration of Eqs. (|2.8|) — (|2.10|) over the spin variables yields the equations and boundary 
conditions for the functions N(x) = JdcN(x,c) and J(x) = JdcJ(x,c): 
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The solution 



AN dJ 

N(0) + | J(0) = 1, iV(L) - | J(L) = 0. (2.12b) 



■J(*)= (| + f) 1 (2-13) 



determines the denominator of Eq. (|2.11 ). 



To determine p we multiply Eqs. ( |2.8| ) and ( |2.9|) with XaX*p an d integrate over c. (Recall 
that xi = c i + «C2, X2 = C3 + ZC4.) It follows upon partial integration that 



y dc XaX%~Q c {M cf) = - ^(SapSfr ~ SapSpa) J dcx P xtf, (2.14a) 

J dc XaXp-^: (Micf) = - ^{T ap 5 Pa - 5 ap Tp a ) J dcxpxtf, (2.14b) 

for arbitrary functions f(x,c). The 2x2 matrices S,T are defined by 

S = itUBT(j z , T — — -j- (a z cos rj — a x sin rj) . (2-15) 

In this way we find that the moments 

N a/3 (x) = J dcXaXp N (%, c), (2.16a) 

Jap(x) = J dcXaXpJfac), (2.16b) 

satisfy the ordinary differential equations 

I , = X! (Tapdpa — dapTpa) N pa + (S a pSpa ~ S ap Sp a ) J pa ~ Jaf3, (2.17a) 
p,cr p,cr 



dx /in ,,.,7 



2 X! {Sapper ~ SapSpa) N pa + ^ {T ap S/3 a — 5 ap Tp a ) J pa , (2.17b) 



with boundary conditions 



N a/3 (x = 0) + | J a/3 (x = 0) = SaxSpt, (2.18a) 
iV a/3 (x = L) - jJ a p{x = L) = 0. (2.18b) 



The mean polarization p is determined by J a p according to 

J(L) ~ V2 



P = Jn(L) "f 22(L) = (? + t ) M - ML)] . (2.19) 



Since Eq. (|2.17Q is linear in the 8 functions N a p (x), J a /3(x) (a, [3 = 1,2), a solution 

requires the eigenvalues and right eigenvectors of the 8x8 matrix of coefficients. These can 
be readily computed numerically for any values oiLjl and uj^t. We have found an analytic 
asymptotic solution for L/£ 3> 1 and u b t ^> (f£/L) 2 , given by 

k , 2-nf sinr? .„ 

p=— — , k= J 1 . 2.20 
smh k i + (2u b t) 2 



In Fig. we compare the numerical solution (solid curve) with Eq. (|2.20|) (dashed curve) for 
L/l — 25 and r\ = 7r/3, / = 1. The two curves are almost indistinguishable, except for the 
smallest values of uj-qt. 




FIG. 2. Average and variance of the spin polarization p of the current transmitted through 
a two-dimensional region of length L = 251, as a function of ujbt, for a magnetic field given by 
Eq. ( |2.l| ) with rj = ir/3 and / = 1. The data points result from Monte Carlo simulations of the 
Boltzmann equation (|2,5|), the solid curves result from the diffusion approximation (p. 7]), and the 
dashed curves are the asymptotic formulas ( 2.20f ) and ( |2.27| ). Notice the transient regime (A), the 
randomized regime (B), and the adiabatic regime (C). 



In a similar way, we compute the second moment of P(p) by multiplying Eqs. (|2.8|) and 



( p.9| ) with XaXpX-yXs an d integrating over c. The result is a closed set of equations 

-N af3y s = {paPyS-N iwpo + ^ap^S JfJkvpa) ~ J a /3~f8, (2.21a) 



£-j-J a {3~{8 — (ZKap^sN^pv + L^gJ^p^ , (2.21b) 



dX I 
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where we have defined 



-^"0/375 — S ol p j 5f3 l/ 5 1 p5s CJ — 5cipSf3 U 5 lp 5s a + 5 Q7i 5 / g ! /S' 7 p5 ( $ CT — 5 a p5p v 5 lp Ss a} (2.22a) 

^affyS = TapSp^SrypSscj — & apX ~i p& + ^afJ,^/3vTy p 5s a — Sa^faSsypTfo, (2.22b) 



N af3l s(x) = J dcXaX*pX-YX*s N ( x i c )> (2.23a) 
J a /3ys(x) = J dcx a XpX~fX*5 J (x,c). (2.23b) 
The boundary conditions on the functions N a ^s an d J a /3-yS are 

7T 

Nafrsix = 0) + -J a /3 7 5(X = 0) = ^i^i^xfe, (2.24) 
7T 

Naf3~/6(X = L) - -Ja/3j6{X = L) = 0. (2.25) 

The second moment p 2 is determined by 

V 1 = (~ + -j) [Juu {x = L)- J im (x = L) — J 22U (x = L) + J 2222 (x = L)] . (2.26) 



.2 

The numerical solution is plotted also in Fig. ||], together with the asymptotic expression 

Var p=i + 2k "f . - J*. (2.27) 

3 3smh(kV3) sinh 2 fc 1 ; 

It is evident from Eqs. (|2.20|) and ( [2.27| ), and from Fig. |2|, that the regime with p = 1, 
Varp = is entered for uj-qt > / [for sin 77 = 0(1)], in agreement with Stern's criterion 
( |1.1| ) for adiabaticity. For smaller uj-qt adiabaticity is lost. There is a transient regime 
^b t ^ {f^/L) 2 , in which the precession frequency is so low that the spin remains in the same 
state during the entire diffusion process. For (f£/L) 2 <C oj-qt <C / the average polarization 
reaches a plateau value close to zero with a finite variance. For a sufficiently non-uniform 
field, / sin 7] 3> 1, we find in this regime p = and Varp = 1/3, which means that the 
spin state is completely randomized. The transient regime, the randomized regime, and the 
adiabatic regime are indicated in Fig. ||] by the letters A, B, and C. 



C. Comparison with Monte Carlo simulations 

In order to check the diffusion approximation we solved the full Boltzmann equation by 
means of a Monte Carlo simulation. A particle is moved from x = over a distance l\ in the 
direction 1; then over a distance £ 2 in the direction 2 , and so on, until it is reflected back 
to x = or transmitted to x = L. The step lengths ti are chosen randomly from a Poisson 
distribution with mean £. The directions 0j are chosen uniformly from [0, 2ir], except for the 
initial direction <fii, which is distributed oc cos0i. The spin components are given by 
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(2.28) 



To find p n , one has to average (|xi| 2 ~~ IX2I 2 )™ over the transmitted particles. The results 
for L/£ = 25 are shown in Fig. |2] (data points). They agree very well with the results of the 
previous subsection, thus confirming the validity of the diffusion approximation for L/£ ^> 1. 

III. WEAK LOCALIZATION 
A. Formulation of the problem 

We turn to the effect of the non-uniform magnetic field on the weak-localization correc- 
tion of a multiply-connected system. We consider the same geometry as in Fig. 0, but now 
with periodic boundary conditions — to model a ring of circumference L. Only the effects of 
the magnetic field on the spin are included, to isolate the Berry phase from the conventional 
Aharonov-Bohm phase. As in the previous subsection, we assume that the orbital motion 
is independent of the spin dynamics. We follow LSG in applying the semiclassical theory of 
Chakravarty and Schmidt! to the problem, however, we start at the level of the Boltzmann 
equation — rather than at the level of the diffusion equation — and make the diffusion 
approximation at a later stage of the calculation. 

The weak-localization correction AG to the conductance is given by 

AG — -~4r I dte- t/T *C{t), (3.1) 
TtnL Jo 

where t v is the phase coherence time and the diffusion coefficient D = v l/d in d dimensions. 
(In our geometry d = 2.) The "return quasi-probability" C(t) is expressed as a sum over 
"Boltzmannian walks" R(t) with R(0) = R(t), 

C(t) = WTr(U + U~). (3.2) 

{R(t)} 

Here VF[R(t)] is the weight of the Boltzmannian walk for a spinless particle. The 2x2 
matrices ?7 ± [R(/f:)] are defined by 

f/ ± = Texp(±^ ['dt'BfKU'^-a}, (3.3) 
I 2/i io J 

where T denotes a time ordering. The factor Tr (U + U~) in Eq. ( ^.2[ ) accounts for the phase 
difference of time- reversed paths. 

The Cooperon can be written in terms of a propagator x, 




that satisfies the kinetic equation 
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+ B j Xafrs( x , Zil 0, 0;; t) - % -^- [( B ( X ) ■ ^ckcAt' ~ <W(B(s) ■ o)^ 



Xa'fS-y'S 



5(t)5(x - X{)5{4> - (pijSa^s- (3.5) 



The Boltzmann operator B is given by 



B = v cos 0— H / — . 

ax r r jo 27r 



The propagator \ is a moment of the probability distribution P(x, <p, U + , U ,t). 

Xaf3~f8 = J 

that satisfies the Boltzmann equation 

with initial condition 

P(x, 0, 17+ LT, 0) = 5(x - Xi)5((f) - (j)i)5{U + - 1)8(U- - 11). 



(3.6) 



(3.7) 



— B — 

dt + + dU+ 



(3. 



(3.9) 



The notation dll + or dU~ indicates the differential of the real and imaginary parts of the 
elements of the 2x2 matrix U + or U~ . We will write this in a more explicit way in the 
next subsection. 

The Boltzmann equation ( |3.8| ) has the same form as that which we studied in Sec. O. 
The difference is that we have four times as many internal degrees of freedom. Instead of a 
single spinor £ we now have two spinor matrices U + and U~ . A first doubling of the number 
of degrees of freedom occurs because we have to follow the evolution of both spin up and 
spin down. A second doubling occurs because we have to follow both the normal and the 
time-reversed evolution. 



B. Diffusion approximation. 

We make the diffusion approximation to the Boltzmann equation ( |3.8| ), by following the 
steps outlined in Sec. ||. The 4x2 matrix u containing the real and imaginary parts of 



u~ 



Re Uti 
Vim 



Re Uf 2 \ 



Re Uf 2 



has a time evolution governed by 



T- 



~d7 



±Z(x)u ± , 



(3.10) 



(3.11a) 
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Z(x) = (jJ-qT 





COS Tj 



cos 77 sin 77 sin — sin 77 cos \ 



sin 77 cos 
sin 77 cos 2^ 

271- f t • • 27r f r. 

\ sin 77 cos — ^— — sin 77 sin — ^— — cos 77 



sin 77 sin ^ L 



2ir fx 

sin 77 sin — j— 

cos 77 




. (3.11b) 



The Boltzmann equation 



becomes, in a more explicit notation, 



d dP d d 

t— P(x, 0, u + , u ,t) = -i cos <P— - E ^+ Z ik^H j P + E -f£P Z ik(x)u kj P 



P + 



dx 

r**d<f/ 
lo 2tt 

We now make the diffusion ansatz in the form 

r27T 



1. 



P(x, (j)' ,u + ,u ,t). 



poo rzn 

/ dte~ t/T * / d(p i P = N + J 
Jo Jo 



COS 1 



(3.12) 



(3.13) 



By integrating the Boltzmann equation over <f), once with weight 1 and once with weight 
cos 0, we obtain two coupled equations for the functions N(x, u + ,u~) and J(x, u + ,u~). Next 
we multiply both equations with U^U~ 5 and integrate over the real and imaginary parts of 
the matrix elements. The moments N a ^s and J a (3~f& defined by 



N aPlS {x) = JdU + J dU- U+ p U~ s N, 
Ja^six) = JdU + J dU- U+ p U; s J : 

are found to obey the ordinary differential equations 
dN a p 1& _ igfiBT 



(3.14a) 
(3.14b) 



dx 



2% 



E [( B ( X ) ' Cr ) QQ '^77' ~ <W(B(x) ■ <T) 77 , J a <fa>6 



o',7' 

— (1 + t/t v ) Jap^S, 



dx 



Ti 



E [( B ( x ) ' ""LcAy - <W(B(x) • <r) 77 , 



-{2T/T v )N a p l8 + 2r5 Q/3 5 75 5(a; - x ; ). 
The periodic boundary conditions are 

^^(O) = N a/3y s(L), Jap-ysi®) = Jaf3-ys(L). 



(3.15a) 
(3.15b) 



(3.16) 

The Cooperon C and the propagator \ of Eqs. (|3.4|) and ( |3~77l) are related to the density iV 
by 



/•oo \ r2TT 

N aPlS {x) = / dte- t/T *— / d, 

JO Z7T JO 

roo 

E^a(^i) = / dte^C(t). 

„n JO 



2tt 



(3.17) 
(3.18) 
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Hence the weak-localization correction ( |3.1|) is obtained from iV by 

AG = -^E Ar ^( a; i)- ( 3 - 19 ) 

The transformation to the local basis of spin states ( |2.3|) takes the form of a unitary 
transformation of the moments N and J, 

N a/ 3 7 s = E Q aa ,Q^yN a ,^y SI Qp lf3 Ql, s , (3.20a) 

a',/3',7',5' 

Jaf3-)S = X/ Qaa'Q-yyJa'P'yS'Q^'flQl'Si (3.20b) 
a> ,/3> ,-y> ,6' 

/ e infx/L cog 2 e -infx/L g j n V \ 

Q(z) = ^ _ e in fx /L gin | g-^/./i cog | J • (3.20c) 
The transformed moments obey 



AN. 



dx 



f3<y'S 



(1 + t/t^) J af3lS , (3.21a) 



- 2 E (S'aa'cVy' — ^ao'S^y) N a 'j3^'S + ^ (T aQ '5 77 / + <5 QCr 'T 77 /) J a 'pys 



d t 

a', 7' a' ,7' 



-(2T/T v )N a p yS + 2r(5 Q/3 5 75 5(x - X;), (3.21b) 

with the same 2x2 matrices S and T as in Sec. [TI|. Because the transformation from N to N is 
unitary, the weak-localization correction is still given by AG = —(e 2 D/iihL) Y, a ,(3N a pp a (x\), 
as in Eq. fl3~T9l) . 

We have solved Eq. ( |3.21|) with periodic boundary conditions by numerically computing 
the eigenvalues and (right) eigenvectors of the 8x8 matrix of coefficients. The resulting 
AG is plotted in Fig. |3] as a function of the tilt angle r\. In the adiabatic regime oj-qt 3> / 
we find the conductance oscillations due to the Berry phase. These are given by@ 

AG = - — ^ ^ML/W) (3 22) 

nh L cosh(L/L ¥ ,) — cos (27r/cos?7) 



analogously to the Aharonov-Bohm oscillations.u (The length = \J Dt^ is the phase- 
coherence length.) In the randomized regime (ft/L) < uj b t < / there are no conductance 
oscillations. Instead we find a reduction of the weak-localization correction, due to dephasing 
by spin scattering. In the transient regime ujbt <C (f£/L) 2 the effect of the field on the spin 
can be ignored,^ and the weak-localization correction remains at its zero-field value 

AG = ~ht co t a nh ( — ] . (3.23) 
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FIG. 3. Weak-localization correction AG of a ring in a spatially rotating magnetic field, as a 
function of the tilt angle r]. Plotted is the result of Eq. ( 3.21 ) for / = 5, L = 500 £, L v = 125 £. 
The upper panel is for uj^t <C 1. From top to bottom: ujbt = 10~ 5 , 10~ 4 , 2 • 10~ 4 , 3 • 10~ 4 , 5 • 10~ 4 , 
10~ 3 , 10 -2 . At ujb t — (f£/L) 2 , the weak-localization correction crosses over from the transient 



regime A of Eq. ( 3.23[ ) to the randomized regime B of Eq. ( |3.30| ). The lower panel is for uj-qt ^ 1. 
From bottom to top: ujb t = 0.1, 1, 2, 5, 10, 100. Here the weak-localization correction reaches the 
adiabatic regime C of Eq. (3.22). 



C. Comparison with Loss, Schoeller, and Goldbart 



If we replace the Boltzmann operator B in Eq. ( |3.5| ) by the diffusion operator —Dd 2 /dx 2 
and integrate over and 4>u we end up with the diffusion equation studied by LSG, 



U = + ^ [B{x) • * x - B(x) • <t 2 ] , 



dx 2 



2h 
1 



■2jt 



2tt 



2tt Jo Jo 

Here cri and cr 2 act, respectively, on the first and third indices of Xa/3-yS- 



(3.24a) 

(3.24b) 
(3.24c) 



The difference between the diffusion equation (|3.24|) and the diffusion equation ( 3.15Q is 
that ( |3.24j ) holds only if wbt <C 1, while ( |3.15| ) holds for any value of ujbt. LSG used Eq. 
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Q3.24j ) to argue that there exists an adiabatic region within the regime uj-qt <C 1. In contrast, 
our analysis of Eq. ( 3.15 ) shows that adiabaticity is not possible if uj-qt <C 1. The argument 
of LSG is based on a mapping of the diffusion equation ( |3.24| ) onto the Schrodinger equation 
studied in Ref. |10[ However, the mapping is not carried out explicitly In this subsection we 
will solve Eq. ( |3.24j ) exactly using this mapping, to demonstrate that the adiabatic regime 
of LSG is in fact the randomized regime B. This mis-identification perhaps occurred because 
both regimes are stationary with respect to the magnetic field strength (cf. Fig. |2|). However, 
Berry-phase oscillations of the conductance are only supported in the adiabatic regime C, 
not in the randomized regime B (cf. Fig. |3|). 

We solve Eq. ( |3.24j ) for the weak-localization correction 



AG 



e 2 D 



a,f3 



n 



(3.25) 



where we introduced the basis of eigenstates \x, a, (3) (with a, (3 = ±1) of the position 
operator x and the spin operators a\ z and a 2z - The operator 7i commutes with 



J 



L d 



2ni dx 2 



+ If {°lz + °2z 



(3.26) 



It is therefore convenient to use as a basis, instead of the eigenstates \x, a, the eigenstates 
\j, a, (3) of J, o"!^, and a 2z . The eigenvalue j of J is an integer because of the periodic 
boundary conditions. The eigenf unctions are given by 



(x,a',f3'\j,a,p) 



exp 



2nix 
L 



U - hfa - \f(3) 



(3.27) 



In the basis {\j, 1, 1), \j, 1, —1), \j, —1, 1), \j, —1, —1)} the operator Ti has matrix elements 

(j',a',{3'\H\j,a,{3) 



-D 



2tt\ 2 

T) ' 



3 3 



( 



\ 





(0 


-f? o 





^ 






J 2 










'3 







f 









\ 





o (j 


+ ff) 









sin 7/ 


— sin 77 







sinr; 


—2 COS 7/ 





— sinr/ 


sin?7 





2 cosr; 


sin r] 







— sin r) 


sin 7] 





/ 



(3.28) 



AG 



Substitution into Eq. (|3.25|) yields 
e 2 D 1 



7fh L 2 

e^_L_ 



E E y><x>P 



a,P 3- 
oo 



j, (3, a 



£ [(7 + f)\f 2 + 7 + j 2 ) + ^ 2 (3/ 2 + 4 7 + 4j 2 + f 2 cos 277) 



j=-oo 
•2\2/ {A 



x [(7 + mr + 2/ 2 7 + 7 2 - 2/ 2 j 2 + 2 7 j 2 + j 4 ) 
+ 2k 2 (/ 4 + 3/ 2 7 + 2 7 2 - / 2 J 2 + 4 7 j 2 + 2/ + / 2 (/ 2 + 7 - 3j 2 ) cos 2 V ) 



(3.29) 
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We abbreviated k = 2ujbt(L/2ti£) 2 and 7 
analytically for k ^> 1, with the result 



(L/2nL ip ) 2 . The sum over j can be done 



AG 



4a_ + 47 + (3 + cos 2r])f 2 4a+ + 47 + (3 + cos 2r/)/ 2 



/a~tan7i\/al 



Q = [/ 4 (g cos 2 2?y - 2 cos 2r]-7)- 32 7 / 2 (l + cos 2r?) 
a± = -7 + i( 1 + 3cos2 ^)/ 2± ^- 



/a^tan7r A 
1/2 



(3.30a) 

(3.30b) 
(3.30c) 



We have checked that our solution ( |3.29|) of Eq. ( p.24|) coincides with the solution of Eq. 
( |3.15| ) in the regime uj-qt <C 1. (The two sets of curves are indistinguishable on the scale of 
Fig. |3|.) In particular, Eq. ( |3.30|) coincides with the curves labeled B in Fig. demonstrating 
that it represents the randomized regime - without Berry phase-oscillations. 



IV. CONCLUSIONS 

In conclusion, we have computed the effect of a non-uniform magnetic field on the spin 
polarization (Sec. [D]) and weak-localization correction (Sec. |III|) in a disordered conductor. 
We have identified three regimes of magnetic field strength: the transient regime uo-qt <C 
(f£/L) 2 , the randomized regime (f£/L) 2 <C uj-qt <C f, and the adiabatic regime uj b t 3> /. 
In the transient regime (labeled A in Figs. [2] and ^), the effect of the magnetic field can be 
neglected. In the randomized regime (labeled B), the depolarization and the suppression 
of the weak-localization correction are maximal. In the adiabatic regime (labeled C), the 
polarization is restored and the weak-localization correction exhibits oscillations due to the 
Berry phase. 

The criterion for adiabaticity is u^tc 3> 1, with u B the spin-precession frequency and t c a 
characteristic timescale of the orbital motion. We find t c = r, in agreement with Stern,! but 
in contradiction with the result t c = t{L/C,) 2 of Loss, Schoeller, and Goldbart.! By solving 
exactly the diffusion equation for the Cooperon derived in Ref. |5], we have demonstrated 
unambiguously that the regime which in that paper was identified as the adiabatic regime, 
is in fact the randomized regime B — without Berry-phase oscillations. 

We have focused on transport properties, such as conductance and spin-resolved trans- 
mission. Thermodynamic properties, such as the persistent current, in a non- uniform mag- 
netic field have been studied by Loss, Goldbart, and Balatsky00 in connection with Berry- 
phase oscillations. These papers assumed ballistic systems. We believe that the adiabaticity 
criterion uj-qt 3> 1 for disordered systems should apply to thermodynamic properties as well 
as transport properties. This strong-field criterion presents a pessimistic outlook for the 
prospect of experiments on the Berry phase in disordered systems. 
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